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Empirical observations have shown that cooperative partners can compete for common 
resources, but what factors determine whether partners cooperate or compete remain unclear. 
Using the reciprocal fig -fig wasp mutualism, we show that nonlinear amplification of interfer- 
ence competition between fig wasps — which limits the fig wasps' ability to use a common 
resource (i.e. female flowers) — keeps the common resource unsaturated, making cooperation 
locally stable. When interference competition was manually prevented, the fitness correlation 
between figs and fig wasps went from positive to negative. This indicates that genetic relatedness 
or reciprocal exchange between cooperative players, which could create spatial heterogeneity or 
self-restraint, was not sufficient to maintain stable cooperation. Moreover, our analysis of field- 
collected data shows that the fitness correlation between cooperative partners varies stochasti- 
cally, and that the mainly positive fitness correlation observed during the warm season shifts to 
a negative correlation during the cold season owing to an increase in the initial oviposition effi- 
ciency of each fig wasp. This implies that the discriminative sanction of less-cooperative wasps 
(i.e. by decreasing the egg deposition efficiency per fig wasp) but reward to cooperative wasps by 
fig, a control of the initial value, will facilitate a stable mutualism. Our finding that asymmetric 
interaction leading to an indeterminate fitness interaction between symbiont (i.e. cooperative 
actors) and host (i.e. recipient) has the potential to explain why conflict has been empirically 
observed in both well-documented intraspecific and interspecific cooperation systems. 

Keywords: asymmetric cooperation; tragedy of the commons; mutualism; 
chaotic oscillation; fig— fig wasp; interference competition 
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1. INTRODUCTION 

Explaining the evolution of cooperation and how coopera- 
tive systems achieve stability remains one of the most 
debated problems in both the biological and the social 
science communities [1,2]. The empirical observations 
and analysis have shown that cooperative actors, in fact, 
take a mixed strategy of either cooperating with the reci- 
pient or cheating/competing with the recipient in both 
intraspecific (e.g. [3-5]) and interspecific cooperative sys- 
tems (e.g. [6-9]). One of the most interesting questions is 
how cooperation can be stably maintained as conflict 
increases between the recipient and the cooperative 
actors under the number increase of involved partners 
or the further action repeat of the cooperative actors in 
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systems with a limited common resource, and the 
increased conflict might lead to breakdown of the coopera- 
tive interaction, even when there is a strong genetic 
relatedness or reciprocal exchange between partners 
[10-13]. Such a breakdown in cooperative behaviour 
represents a 'tragedy of the commons' in the evolution 
of cooperation [14]. 

Classical theories of evolution of cooperation based 
on an assumption of symmetrical interactions argued 
that cooperation could be maintained either through 
the self-restraint of cooperative actors or through 
spatial heterogeneity created by the recipient because 
of high genetic similarity or reciprocal exchange 
between the cooperative partners [10,15-18]. However, 
theoretical explorations showed that taking a non- 
cooperative strategy or even being a parasite within 
these systems will be more advantageous than the coop- 
erating strategy of self-restraint (i.e. not competing with 
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partners) ; after the utilization of the common resource is 
saturated, individuals employing the pure cooperation 
strategy (not competing with others) may not receive 
any fitness benefit [10,19,20]. The spatial heterogeneity 
created by the recipient is also not a credible mechanism 
for maintaining the stability of cooperative systems. 
Recent theoretical advances have shown that heterogen- 
eity will often emerge but fails to promote the evolution 
of cooperation, because the mutant that can overcome 
the spatial heterogeneity to use the common resource 
at the expense of recipient will have higher fitness than 
cooperating actors in a limited local space [21]. 

Empirical observation and data analysis have fur- 
ther shown that self-restraint of cooperative actors or 
spatial heterogeneity created by recipients cannot suffi- 
ciently maintain the stable cooperative interaction in 
both intraspecific and interspecific cooperation systems. 
In almost all of the well-documented intraspecific 
cooperation systems, the literature shows that coopera- 
tive actors (i.e. the subordinates) are in fact competing 
with the recipient (i.e. king or queen) for the common 
resource or tend to be lazy to reproduce by themselves 
(e.g. in lion, see [3]; in mole rodent, see [5]; in bee and 
ant, reviewed by [22]). Explorations of interspecific 
cooperation systems also show that the 'honest' cooperative 
actors of a mutualism (i.e. the symbiont) can compete with 
their host (i.e. cooperative recipient) in several well-known 
mutualism systems (e.g. in fig-fig wasp mutualism, see 
[6,23,24]; in mutualism between moth Greya politella and 
plant Lithophragma parviflorum, see [25]). Our previous 
manipulative experiments with the fig -fig wasp mutual- 
ism have shown that the cooperative interaction between 
the fig and its pollinator wasps will transform into a conflict 
interaction when the unused common resource (unused 
female flowers of figs) decreases [9] . 

That both cooperation and conflict exist in these 
cooperative systems show that pure cooperation equili- 
brium might not exist in cooperative systems, and that 
the systems that display cooperation might be in a non- 
equilibrium state (e.g. a chaotic system, see [26-28]). 
In non-equilibrium systems, the behavioural impact of 
one partner on the other partner is indeterminate, and 
therefore the interaction between involved partners 
cannot be expected to be constant. However, such systems 
with indeterminate behavioural interactions between 
particles/elements/partners could be maintained at a 
steady state through asymmetric interaction (e.g. [29] 
nuclear physics; see [30] thermodynamics; see [31] 
economics). Theoretically, the control of initial pertur- 
bation can favour the steady state in non-equilibrium 
systems [32,33] . In cooperative systems, asymmetric inter- 
action between the cooperative actors and recipients (e.g. 
density-dependent interference competition among coop- 
erative actors, but not among recipients; the recipient 
might discriminatively sanction against non-cooperative 
actors but reward the cooperative actors) might lead to 
a non-equilibrium state, but the control of initial 
common resource utilization efficiency will favour the 
steady-state maintenance of the cooperation systems 
(see mathematical simulation of [20]). 

In this paper, we use direct fitness correlation analysis 
to examine whether the well-documented interspecific 
mutualism between figs and fig wasps is stabilized by an 



asymmetric interaction. In the fig- fig wasp mutualism, 
figs provide pollinator wasps with part of their female 
flowers for oviposition and the eggs of pollinators develop 
into adults in the ovaries of flowers (galls); in exchange, 
wasp pollinators carry pollen from the mature fruits to 
receptive fruits. In the process of oviposition, pollen is dis- 
persed to female flowers by foundresses (pollinators in 
fruit cavities) , and part of fertilized female flowers devel- 
ops viable seeds. In monoecious figs, because both 
pollinator offspring and viable seeds of figs use the same 
resource (i.e. female flowers), there should be conflict in 
using female flowers, which could disrupt the cooperation 
system [6,9,23]. 

It is generally believed that spatial heterogeneity (i.e. 
the structural barrier of flowers) created by host plants 
and/or an evolutionary constraint on fig wasps prevent 
pollinators (symbionts) from over-using common resource 
at the expense of host plants, thus keeping mutualistic 
cooperation in a stable equilibrium state. However, further 
quantitative examination has shown that neither the 
hypothesized spatial heterogeneity (structural barrier of 
flowers such as the heterogeneity in style or pedicel 
length) created by plants nor the evolutionary cons- 
traint of symbionts in the fig- fig wasp mutualism to be 
sufficient to prevent fig wasps from over-using female 
flowers at the expense of viable seeds [9,24,34,35]. 
This is in spite of many suggested mechanisms of spatial 
heterogeneity or evolutionary constraint for symbionts. 

In a cooperative system, if spatial heterogeneity or 
self-restraint can prevent a direct conflict between recipi- 
ent and cooperative actors, only a positive fitness 
correlation (i.e. viable seed-wasp offspring correlation) 
can be observed and a negative fitness correlation (con- 
flict) will not be present even when the use of the 
common resource is saturated [9,10,17]. However, posi- 
tive fitness correlation (i.e. a cooperative interaction) 
can also be observed if the common resource is unsatu- 
rated, given that genetic relatedness or reciprocal 
exchange between the recipient and cooperative actors 
is higher than the cost/benefit ratio [11,12]. In the 
research described in this paper, we used fig (Ficus race- 
mosa) and its obligate pollinator Ceratosolen fusciceps to 
examine how the steady state in a non-equilibrium state 
between figs and their pollinator wasps is maintained 
under the variation of environmental or ecological factors 
through interference competition among foundresses (i.e. 
pollinator wasps in fruit cavities). The interference com- 
petition among the foundresses that could decrease the 
oviposition of the flowers [36], however, might also 
decrease the pollination of the flowers. Following from 
our previous paper [36], in this paper we would like 
to explore the remaining problem on why both 
cooperative (i.e. positive fitness correlation) and conflict 
interactions (i.e. negative fitness correlation) can be 
observed between figs and their pollinator wasps under 
the interference competition among the fig wasps. 

2. MATERIAL AND METHODS 
2.1. Study species 

Ficus racemosa Linn. (Ficus sycomorus) is distributed 
from India to Australia. Ficus racemosa is a large tree 
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that can reach 30 m high and bears cauliflorous fruits 
synchronously within the tree in very large numbers. 
It grows mainly in moist valleys or along rivers. Trees 
of F. racemosa usually grow in groups of five to 10. 
Ficus racemosa is pollinated by the wasp species C. fus- 
ciceps Mayr (Agaonidae). The foundress number per 
receptive fruit of F. racemosa usually range from 1 to 
30, but can sometimes reach more than 70. The aver- 
aged foundress number per receptive syconium varies 
greatly with different habitat sites and seasons [37]. 

2.2. Introduction method 

Interference competition, which should exist if foun- 
dresses (cooperative actors) simultaneously enter fruit 
cavities, can be excluded partly or totally if foundresses 
sequentially enter fruit cavities at time intervals [36,38]. 
We designed an experiment to examine whether there 
was interference competition among pollinator fig 
wasps, and to observe whether there was a negative fit- 
ness correlation when interference competition was 
excluded. In one set of treatments, we simultaneously 
introduced pollinators (manually collected wasps) into 
receptive fruits within 30 min. In these experiments of 
simultaneous introduction of foundresses, we conducted 
six treatments in which two, five, seven, nine, 15 and 20 
foundresses in each treatment were simultaneously 
introduced into each receptive syconium. 

With sequential introduction of foundresses, we con- 
ducted four treatments in which two, five, seven and 
nine foundresses in each treatment were separately 
introduced. In the two-foundress treatment, we intro- 
duced one on the first day at 09.00 h, and the second 
on the second day at 09.00 h for each receptive syco- 
nium. In the other three treatments with, respectively, 
five, seven and nine foundresses in total for each 
receptive syconium, we introduced two on the first 
day at 09.00 h and then sequentially introduced two, 
or one if only one was left, every 3-4 h in the day 
time or 10-12 h in the evening. Foundresses were col- 
lected using an insect net from the air surrounding 
receptive syconia. 

In the above comparison experiments between 
the simultaneous and sequential introductions, the 
data of oviposition of flowers with two, five, seven and 
nine foundresses have appeared in our previous paper 
[36]. Because the interference competition occurs in 
both oviposition and pollination of the foundresses, in 
this paper we would like to see whether there is a differ- 
ence in the intensity of interference competition among 
the foundresses between the oviposition and pollination 
of the flowers. 

Because the oviposition efficiency of each foundress 
(i.e. initial value of oviposition) differs between morning 
and afternoon (table 1), we conducted an extra four 
treatments with both sequential and simultaneous intro- 
duction with nine foundresses to compare whether the 
oviposition efficiency of each foundress affects final fitness 
correlation in the presence of density-dependent interfer- 
ence competition among the foundresses. In these four 
extra experiments, the first introduction was at 06.00 h 
and at 15.00 h, while the others were treated in the 
manner identical to the other introductions. These extra 



four treatments were only conducted in one sample site 
(crop field) with the same two trees to ensure that other 
environmental conditions were the same. 

In the experiment designed to determine oviposition 
efficiency, we manually introduced two foundresses 
into receptive fruits. Syconia pollinated by a single 
foundress often abort [40] , and there is limited interfer- 
ence competition for pollination or oviposition between 
foundresses in large figs (e.g. those of F. racemosa) 
when foundress number is two [38]. The mean pollina- 
tion and oviposition per foundress for syconia- 
containing two foundresses can therefore be used to 
evaluate oviposition efficiency of foundresses. We con- 
ducted such experiments in both the warm and cold 
seasons, and morning and afternoon within the 
same day. 

In these experiments, all of the treated fig fruits were 
enclosed in nylon bags during development to prevent 
any oviposition by other wasps. In the above exper- 
iments, we selected two trees in a crop field and 
another two trees in a locally fragmented forest. 

2.3. Natural data collection 

We also counted the number of viable seeds in each fig, 
the number of fig wasp offspring (galled flower number) 
and the amount of unused common resource (vacant 
female flowers) over years from four sampling sites at 
the centre of the primary forest, edge of the primary 
forest, within a locally fragmented forest and from iso- 
lated fig trees by rubber trees, roads or crop fields. We 
collected these data over most months of the year. Ovi- 
posited flowers (galls), viable seeds and vacant female 
flowers were counted in premature syconia (pre-D 
phase) or mature syconia (D phase). We vertically cut 
each syconium into eight slices passing through the 
ostiole, and then randomly selected two or three slices 
to count all galls, seeds and vacant flowers. Owing to 
the difficulty of counting all vacant flowers, only galls 
and seeds of the remaining five or six slices were 
counted. We then calculated the percentage of the 
developed flowers (galls + seeds) per syconium and esti- 
mated the total number of the female flowers per 
syconium using the following calculation: the total 
number of the female flowers = (total galls + total 
seeds) /proportion of the developed flowers. 

2.4- Correlation analysis 

Cooperation or conflict interaction can essentially be 
described by the direct or inclusive fitness correlation 
coefficient between involved players in both intraspecies 
and interspecies cooperation systems [1,9,20,41,42]. 
However, if the cooperation interaction of involved 
players is not in a stable equilibrium (e.g. chaotic oscil- 
lation) a constant fitness relationship between involved 
players cannot be expected, possibly presenting an inde- 
terminate interaction between the individual players. 
On the other hand, the whole system could still be in 
a stable state (e.g. [32,33,43]). In the case where there 
is no constant relationship between interacting part- 
ners, the fitness regression coefficient using simple 
linear regression or any other parametric correlation 
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Table 1. Factors that might influence oviposition efficiency of foundresses (oviposition efficiency: galled flowers per foundress), 
(mean + s.e.). The experimental data showed that the lifespan of the pollinator was strongly affected by the temperature and 
humidity of environmental conditions [24,39], which might result in the oviposition efficiency difference between morning and 
afternoon, or between warm and cold season. The data in this table also show that the fruit diameter will significantly affect 
the oviposition efficiency of each foundress, conforming to the correlation analysis result in §3. Here, O/F is the offspring (gall) 
number per foundress, and FD is the fruit diameter (millimetre). 



factors 



between different fruit 
diameters (nine 
foundresses) 



between different fruit 
diameters (two 
foundresses) 



between morning 
and afternoon (two 
foundresses) 



between warm 
and cold seasons 
(two foundresses) 



treatment 1 



treatment 2 



comparison 
between two 
treatments 



O/F: 114.1 + 5.4 
(n=20) 

FD: 41.2 + 1.1 

O/F: 155.9 + 12.2 
(n=20) 

FD: 50.3 + 1.0 

O/F: ^ = 3.1; d.f. = 38, 
p<0.01;FD: *=6.3, 
d.f. = 38, p < 0.001 



O/F: 264.5 + 11.2 
(n=21) 

FD: 45.8 + 0.6 

O/F: 273.4 + 9.0 
(n=20) 

FD: 36.4 + 0.2 

O/F: £=0.6; d.f. = 39, 
p > 0.05; FD: t= 14.9, 
d.f. = 39, p < 0.001 



O/F: 502.3 + 13.6 
(09.00 h, n=20) 



O/F: 229.9 + 9.0 
(15.00 h, n=20) 

t= 16.7, d.f. = 38, 
p < 0.001 



O/F: 381.8 + 30.4 
(n= 34; Apr., Mar.. 
Sep.) 

O/F: 812.1 + 56.9 
(n= 34; Nov., Dec.) 

^ = 6.7, d.f. = 66; 
p < 0.001 



Table 2. The correlation coefficient between viable seeds and the galled female flowers in different conditions when the 
foundress number is 9, controlling for the syconium size and the total female flower number. The results here showed, given the 
foundress number, that the correlation coefficient between viable seeds and galled female flowers varies greatly under different 
introduction treatments, as well as different environmental conditions. 



sample site and 


sample 


introduction 


total female flowers per 


syconium size 


correlation 


introduction time 


size 


treatment 


syconium (Mean + s.e.) 


(Mean + s.e.) 


coefficient 


samples from different habitat sites and different introduction treatments 






crop field (at 09.00 h) 


20 


simultaneous 


5443.1 ± 161.8 


50.3 ± 1.0 


0.74*** 


crop field (at 21.00 h) 


21 


sequential 


5815.9 ± 295.8 


50.4 + 0.8 


-0.67*** 


fragmented forest (at 


20 


simultaneous 


4084.6 + 81.1 


41.2 + 1.1 


0.68** 


21.00 h) 












fragmented forest (at 


20 


sequential 


4266.1 + 82.7 


38.6 + 0.2 


-0.07 n.s. 


21.00 h) 












samples from the same 


habitat site with different treatments 






crop field (at 06.00 h) 


20 


simultaneous 


4886.2 + 96.9 


46.1 + 0.3 


0.06 n.s. 


crop field (at 06.00 h) 


20 


sequential 


4791.8 + 67.3 


46.0 + 0.3 


-0.96*** 


crop field (at 15.00 h) 


20 


simultaneous 


4701.5 + 98.0 


45.5 + 0.3 


0.75** 


crop field (at 15.00 h) 


20 


sequential 


4966.8 + 91.0 


46.9 + 0.4 


-0.64** 



n.s. is not significant at p = 0.05. 

**p<0.01. 

*** p< 0.001. 



analysis will be the averaged value of their oscillating 
correlation coefficients, and therefore would not 
describe the true relationship between the involved 
players [9,43]. In an indeterminate system (e.g. chaos 
system) , the direct fitness or inclusive fitness correlation 
coefficient (i.e. fitness regression coefficient) between 
recipient and cooperative actors should vary or oscillate 
with variation of affecting factors, which could be the 
frequency of cooperative actors, or common resource 
availability or other factors [9,20,27,43]. 

In the manipulated experiments, because the effect 
of indeterminate environmental or ecological factors 
was mostly prevented, we used a simple linear 
regression to analyse the correlation between the 
viable seeds and egg deposition of the foundresses (i.e. 
galls), when the foundresses were simultaneously 



introduced into the syconium cavity. When the foun- 
dresses were sequentially introduced into the 
syconium cavity, we used quadratic function to fit the 
distribution pattern of the viable seeds as a function 
of wasp offspring. This analysis was able to determine 
whether any positive contribution of egg deposition 
translated to a negative effect on seed production. 

Because the correlation coefficient between the 
viable seeds and wasp offspring number varies so 
greatly a constant correlation coefficient cannot be 
expected under the variation of the environmental 
change, even in the experiments given the foundress 
number (table 2). Under natural conditions, the factors 
that will affect the fitness correlation between fig and fig 
wasps will be indeterminate, and therefore the linear 
regression or generalized linear regression which 
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assumes that the function response between the interacted 
variables could be given cannot be used to estimate such 
indeterminate interaction between fig and fig wasps. For 
the data collected from the natural condition, we used 
varying correlation coefficient with non-parametric esti- 
mation to describe the correlation between viable seeds 
of figs and wasp offspring of fig wasps as a function of 
other factors. Through such a method, we can describe 
how the fitness interaction between figs and fig wasps 
varied with changes in other factors. 

The varying coefficient method with non-parametric 
estimation does not pre-specify functional forms of the 
involved variables, but uses information purely from 
data to estimate a curve/function. Such varying coeffi- 
cient method is obviously more flexible and robust to fit 
data than any parametric method, because if a para- 
metric structure of the model is pre-assumed, there is 
a risk of the wrong model and conclusions [44,45]. 
When the functional response between the interacted 
variables can be given to be a specific form, the varying 
coefficient is then reduced to the traditional parametric 
correlation analyses, such as linear regression or gener- 
alized linear regression methods. As such, the varying 
coefficient is more general and more flexible for use in 
empirical data analyses, and then, different from the 
traditional methods, it might help us identify whether 
or not species interaction varies with respect to other 
factors, and what factors will be crucial to maintain 
or change the species interaction. For details of the 
varying coefficient analysis model and the comparison 
between the varying coefficient method and traditional 
linear regression or generalized linear regression see the 
electronic supplementary material, online appendix or 
the paper by Shi et al. [46]. 



3. ANALYSES AND RESULTS 

Similar to the results of Bronstein et al. [38] , there was no 
interference competition when there were two foun- 
dresses, and both viable seed production and wasp 
offspring (gall) production were not significantly differ- 
ent between treatments of simultaneous and sequential 
entry of fig wasps (figure la, 6). In contrast, there was 
lower oviposition and pollination with simultaneous 
entry of pollinators in treatments with five, seven and 
nine foundresses than with sequential entry. The mean 
comparison with general linear model including the 
syconium size and the total number of female flowers as 
covariates for seed production were for five foundresses, 
F 1 78 = 11.116, p < 0.001; for seven foundresses, F\ 79 — 
37.912, p< 0.001; and for nine, F lj79 = 47.979, p < 
0.001. For gall production, results were for five foun- 
dresses, Fi ? 7 8 = 44.332, p < 0.001; for seven foundresses, 
F 179 = 66.369, p < 0.001; and for nine, F^ 79 = 67.910, 
p < 0.001. The intensity of interference competition on 
the production of galled flowers is higher than on the 
viable seeds (figure lc). 

With simultaneous entry of pollinators, the nonlinear 
amplification of interference competition among the 
foundresses cancelled out extra pollination and oviposi- 
tion after the number of foundresses reached a threshold 
(figure 2). The interference competition among 
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Figure 1 . The effect of foundress abundance on the intensity of 
interference competition between the foundresses on the gall 
and seed production in F. racemosa. (a) Gall production 
(mean + s.d.), which is similar to fig. 5 of Wang et al. 
[36], (b) viable seed production (mean + s.d.), (c) interference 
competition intensity ((number of sequential introduction — 
number of simultaneous introduction) /number of sequential 
introduction) . That the efficiency of pollination per foundress 
is higher than oviposition (i.e. initial value) is selected by the 
discriminative sanction of less cooperative or cheating wasps 
but reward to the cooperative wasps by host ([40]; R.-W. 
Wang et al. 2009, unpublished data), (a, b) Light grey bars, 
foundress sequential entry; dark grey bars, foundress simul- 
taneous entry, (c) Diamonds with solid line, galls; squares 
with solid line, seeds. 
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Figure 2. The production of wasp offspring and viable seeds as 
a function of the number of simultaneously entering foun- 
dresses (mean + s.e.). The data were collected from two 
trees in a crop field and sample size of each treatment greater 
than 20. Diamonds with solid line, seeds; squares with solid 
line, galls. 
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Figure 3. Proportion of seeds (seed number/total female flower number) as a function of the proportion of wasp offspring (galled 
flower number/total female flower number) in F. racemosa with different treatments, (a, b) Collected from crop field, (a) Foun- 
dresses introduced simultaneously; the foundress numbers were two, three, five, seven and nine with sample size 21, 20, 20, 20 and 
20, respectively (linear: n = 101, r= 0.67, p < 0.001). (b) Foundresses sequentially introduced; foundress numbers were one, two, 
five, seven and nine with sample size 21, 20, 20, 23 and 65, respectively (quadratic: Ftest = 201.87, p < 0.001). (c,d) Collected 
from locally fragmented forest, (c) Foundresses introduced simultaneously; foundress numbers were two, three, five, seven and 
nine with sample size 20 for each treatment (linear: n= 100, r=0.52, p< 0.001). (d) Foundresses sequentially introduced; 
foundress numbers were one, two, three, five, seven and nine with sample size 20 samples for each treatment (quadratic: 
Ftest = 23.14, p< 0.001). 



foundresses decreased the common resource utilization 
(i.e. developed female flowers) in both pollination and ovi- 
position. Such a decreased utilization of common resource 
from interference competition among wasps makes the 
common resource utilization unsaturated. However, in 
experiments with sequential entry of pollinators, the 
viable seed production was obviously an inversed U- 
shape as a function of wasp offspring production, which 
differed from positive correlation for simultaneous entry 
(figure 3). The correlation coefficient between viable 
seeds and wasp offspring varies greatly in the experiments 
with the change of foundress numbers, introduction 
methods and habitat sites (tables 2 and 3) . In the natu- 
rally collected data, the impact of wasp offspring on 
production of viable seeds (indicated by correlation coef- 
ficients) oscillated with the change in density of wasp 
offspring or the availability of unused common resource 
(vacant female flowers; figure 4). 

The distribution patterns of viable seed production 
as a function of galled flower number differ between 
different sample sites (figure 3). This is expected 
because the interference competition impacts depend 
significantly on fruit diameters and the total number 
of female flowers, and these parameters differ between 
these two sample sites. In the simultaneous introduc- 
tion experiment with nine foundresses, the correlation 
analysis showed that the fruit diameter was positively 
correlated with oviposition efficiency of foundresses 



(galled flower number per foundress), (n = 40, r = 0.47, 
p < 0.05). The oviposition efficiency of foundresses 
were also significantly correlated with the total number 
of female flowers (n=40, r=0.53, p < 0.001). The 
fruit diameter size and the total female flowers of samples 
from the crop field are on average larger than from the 
locally fragmented forest (diameter: t= 25.986, d.f. = 
365, p< 0.001; female flowers: t= 11.175, d.f. = 365, 
p< 0.001). 

In the fig-fig wasp mutualism, the oviposition effi- 
ciency per foundress varies greatly across different 
ecological or environmental conditions. In our exper- 
iments, the foundresses had higher oviposition efficiency 
in the morning than in the afternoon, and in particular 
demonstrated higher efficiency in the cold than in the 
warm season (table 1). This is because foundresses have 
much longer lifespans at lower temperature [24] and 
higher humidity [39]. In the morning, the humidity is 
higher than in the afternoon and the temperature is 
usually lower than in the afternoon. In table 2, the foun- 
dress number is the same in all of the treatments, but the 
fitness interaction between viable seeds and galled flowers 
varies greatly under different foundress introduction treat- 
ments (i.e. simultaneous or sequential introduction), in 
different sample sites with different fruit size or total 
female flower number, and under different initial values 
of female flower utilization by foundresses at different 
introduction times (i.e. morning or afternoon; table 1). 
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Table 3. The correlation coefficient between viable seeds and the galled female flowers in different conditions and foundress 
number, controlling for the syconium size and the total female flower number. The results here showed that the correlation 
coefficient between viable seeds and galled female flowers varies greatly under different introduction treatments and different 
foundress numbers, as well as different environmental conditions. 



sample site and 


foundress 


sample 


introduction 


total female flowers per 


syconium size 


correlation 


introduction time 


number 


size 


treatment 


syconium (mean + s.e.) 


(mean + s.e.) 


coefficient 


crop field 


2 


21 


simultaneous 


4483.0 + 104.5 


45.8 + 0.6 


0.15 n.s. 


crop field (at 


2 


20 


sequential 


4360.7 + 107.6 


44.4 + 0.5 


0.26 n.s. 


21.00 h) 














fragmented forest 


2 


20 


simultaneous 


4136.5 + 66.9 


36.4 + 0.2 


-0.29 n.s. 


/ i o 1 r\r\ 1 \ 

(at 21.00 h) 














fragmented forest 


2 


20 


sequential 


4055.6 + 45.2 


35.9 + 0.5 


0.50 n.s. 


(at 21.00 h) 














crop field (at 


5 


20 


simultaneous 


4698.7 + 86.4 


48.2 + 1.1 


0.29 n.s. 


09.00 h) 














crop field (at 


5 


20 


sequential 


4753.9 + 61.1 


49.5 + 40.9 


-0.03 n.s. 


21.00 h) 














fragmented forest 


5 


20 


simultaneous 


4169.3 + 107.8 


40.2 + 0.6 


0.07 n.s. 


(at 21.00 h) 














fragmented forest 


5 


20 


sequential 


4217.8 + 36.9 


39.3 + 0.8 


-0.20 n.s. 


(at 21.00 h) 














crop field (at 


7 


20 


simultaneous 


4393.8 + 50.7 


50.7 + 0.5 


0.40 n.s. 


09.00 h) 














crop field (at 


7 


23 


sequential 


5578.1 + 260.8 


52.1 + 0.8 


-0.401 n.s. 


21.00 h) 














fragmented forest 


7 


20 


simultaneous 


4075.1 + 68.2 


43.0 ± 0.6 


-0.26 n.s. 


(at 21.00 h) 














fragmented forest 


7 


20 


sequential 


4309.3 + 29.1 


36.0 ± 0.3 


-0.79*** 


(at 21.00 h) 















n.s. is not significant at p = 0.05. 
***p< 0.001. 



The pollination efficiency of foundresses is higher 
than oviposition efficiency of foundresses (figure 1), 
leading to faster rates of pollination saturation than ovi- 
position saturation in the female flowers of figs. When 
female flower utilization is saturated (i.e. flowers are 
either pollinated or oviposited), an increase in egg depo- 
sition leads to a decrease in seed production, and 
therefore a negative fitness interaction could be 
observed. In the warm season, the fitness interaction 
between figs and fig wasps was mainly positively corre- 
lated, but mainly negative in the cold season (figure 4). 
Such seasonally dependent seed -wasp offspring inter- 
action patterns result because the initial rate of egg 
deposition (i.e. egg deposition efficiency of each foun- 
dress) in the cold season is much higher than in the 
warm season (table 1). 



4. DISCUSSION 

The results of our investigation of interspecific 
cooperation between figs and fig wasps show that a 
local cooperative equilibrium between recipient (e.g. 
hosted figs) and cooperative actors (e.g. pollinator 
wasps) can be maintained by asymmetric interference 
competition (i.e. interference competition was mainly 
among the cooperative actors, not directly for the reci- 
pient). The interference competition among the 
cooperative actors decreases the common resource 



utilization by cooperative actors and kept the 
common resource utilization unsaturated, thus main- 
taining the cooperative interaction (i.e. positive fitness 
interaction; figures 2 and 3). However, the impact of 
interference competition among the cooperative actors 
will be greatly affected by ecological or environmental 
constraints (e.g. dispersal barrier, number of coopera- 
tive actors, availability of common resource and 
common resource utilization efficiency of foundresses), 
and such ecological or environmental constraints 
might be indeterminate. The fitness interaction shown 
in figure 4 therefore becomes a stochastic process, and 
the probability of a positive fitness interaction differs 
under different ecological or environmental conditions. 

Our results imply that the common resource of female 
flowers for recipient (host plant) and cooperative actors 
(symbionts) was global rather than local, and direct fit- 
ness conflict between fig plants and pollinator wasps 
would occur after the common resource utilization was 
saturated in some specific situations. For example, in 
the environment with higher humidity or lower tempe- 
rature, the oviposition efficiency can be increased, or 
when the foundresses could enter syconium with a 
longer interval, the interference competition among the 
foundresses can be lessened. This observation suggests 
that genetic relatedness or reciprocal exchange may be 
important (by creating spatial heterogeneity or exer- 
cising recognition discrimination against additional 
cooperative actors or cheating individuals) but is not 
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Figure 4. The Pearson correlation coefficient between the fitness of figs (viable seeds) and the fitness of pollinator wasps (offspring 
number) of cooperating species using a non-parametric estimation method, (a) Samples collected in the warm season (n= 218, 
March- October), (b) Samples collected in the cold season (n= 134, November -February), (c) Fitness correlation coefficient 
between seeds and wasp offspring as function of unused commons (vacant female flowers; n= 352, over years), (d) Fitness cor- 
relation coefficient between seeds and wasp offspring as function of wasp offspring (n = 352, over years). The bandwidth h is 400, 
which was obtained by cross-validation method. 



by itself a sufficient explanation for how cooperation is 
maintained. This is the first direct empirical evidence 
in support of the argument that the Nash equilibrium 
or evolutionary stable strategy in cooperative systems 
may be inaccessible [19,47] and that an explanation of 
cooperation based on spatial chaotic oscillations might 
be more credible [26,27]. The genetic relatedness or reci- 
procal exchange might just be a pivot rather than an aim 
of evolution towards cooperation [28], and the higher 
genetic relatedness or reciprocal exchange will reduce 
the probability of competition [20] . 

Genetic relatedness or reciprocity exchange results 
from limited dispersal from the original colony [17,48]. 
However, limited dispersal might also result in compe- 
tition between involved partners, not only between 
cooperative actors and the recipient, but also among 
the cooperative actors (i.e. pollinator wasps in fig- 
fig wasp mutualism, see [36]) [12]. As is shown in 
figures 1 and 3, more intensive competition between 
the cooperative actors means there will be a lower prob- 
ability of the common resource being over-exploited, 
and therefore a less probability of conflict in the inter- 
action between cooperative actors and the recipient. 
Theoretically, higher dispersal difficulty for cooperative 
actors (creating higher genetic relatedness or reciprocity 



between the recipient and the cooperative actors) 
will increase the interference competition among the 
cooperative actors, and therefore higher rates of co- 
operation will be expected in systems with higher 
dispersal difficulty [20]. The higher dispersal difficulty 
for cooperative actors will make the sanction or repres- 
sion of less-cooperative wasps by host or dominant 
individuals more credible [40]. Through such a control 
of initial value (e.g. by decreasing the egg deposition 
efficiency per fig wasp in fig -fig wasp mutualism), the 
steady state in non-equilibrium of cooperation systems 
might be maintained [32,33]. 

Essentially, the cooperative interaction is maintained 
by density-dependent interference competition (i.e. a 
process of nonlinear amplification) owing to suppression 
of resource consumption, which prevents resource use 
from being saturated (figures 3 and 4c). Variation in 
the availability of a common resource, the intensity of 
interference competition or common resource utiliza- 
tion efficiency will greatly affect the saturation of the 
common resource utilization, whereas the availability 
of unused common resource will directly determine the 
fitness correlation relationship between the symbiont 
and the host (figure 4c). Further, the above variables 
are also greatly affected by ecological or environmental 
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factors. As such, the relationship between the symbiont 
and the host, therefore, might be indeterminate. Our 
discovery of the factors influencing this oscillating 
relationship reinforces previous findings that the inter- 
action between symbionts and hosts might oscillate 
between mutualism and antagonism across temporal 
or spatial variability [24,25,41,49,50]. The interaction 
between the reciprocal mutualists might be conditional 
[51] . This also explains why increased density of coopera- 
tive actors leads to direct conflict (competition) between 
cooperative actors and recipient for a limited common 
resource [9,12]. 

Density-dependent (i.e. nonlinear) interference com- 
petition among cooperative actors might be a general 
mechanism for maintaining cooperation rather than a 
special mechanism of the fig -fig wasp mutualism. 
In intraspecific cooperative systems, there exist dispersal 
barriers for involved partners and a larger number of coop- 
erative actors than recipients [1,3,17], a condition similar 
to that of the fig-fig wasp interspecific cooperation 
[9,48]. Thus, competition among the cooperative actors 
might be more intense than between recipient and coop- 
erative actors, owing to asymmetric interaction between 
recipient and cooperative actors. Such density-dependent 
interference competition among cooperative actors will 
lead to an indeterminate fitness interaction between 
cooperative actors and recipient [20,52]. Different fitness 
interactions might lead to different evolutionary and be- 
havioural strategies. This could explain why cooperative 
actors take mixed strategies of cooperation or competition, 
rather than pure cooperation with the recipient in 
genotype, phenotype or behaviour [3,7,9,53]. 

A remaining question is what dynamics drive the 
players involved to pay evolutionary costs to punish 
cheating or less- cooperative individuals with the poten- 
tial to disrupt the cooperative interaction [12,40,54]. If 
a high probability of cooperative interaction can be 
maintained by environmental or ecological restraint 
through the mechanism of mutual interference among 
cooperative actors, the net interest of the recipient will 
be greater for cooperation than conflict over an entire 
life cycle. It would be advantageous to the recipient 
(i.e. host of mutualism or dominant of intraspecific 
cooperation systems) to pay costs to maintain a stable 
cooperative interaction by punishing less cooperative or 
cheating behaviours of cooperative actors [52,55,56]. 
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